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Abstract 

The application of the hyperspherical harmonic approach to the case of the N — d scattering 
problem below deuteron breakup threshold is described. The nuclear Hamiltonian includes two- 
and three- nucleon interactions, in particular the Argonne vis, the N3LO-Idaho, and the Viou^k 
two-nucleon, and the Urbana IX and N2L0 three-nucleon interactions. Some of these models 
are local, some are non-local. Also electromagnetic effects are included. Accurate calculations 
for many scattering observables at various center-of-mass energies are performed and the results 
are compared with the available experimental data. Furthermore, a analysis of some of the 
Hamiltonian models has been performed to compare their capability to describe the scattering 
process. 
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I. INTRODUCTION 



One of the main inputs for any study on nuclear systems within a non-relativistic frame- 
work is the model used to describe the nuclear interaction, i.e. the nuclear Hamiltonian. 
Nowadays, it is common practice to use, at least for few-nucleon systems, Hamiltonian mod- 
els composed of a two-nucleon plus, for A > 3, a three-nucleon interaction (TNI). The mod- 
ern two-nucleon interaction models have a large number of parameters and can reproduce 
the deuteron properties and the nucleon-nucleon scattering data up to the pion threshold 
with a x^/datum ~ 1. Among them, the Argonne fig (AV18) [1] and the charge- dependent 
Bonn (CDBonn) ji2i] explicitly include charge-symmetry-breaking terms in the nuclear in- 
teraction, in order to reproduce equally well the np and pp data. Recently, a number of 
two-nucleon interaction models have been derived by many authors within an effective field 
theory (EFT) approach, up to next-to-next-to-next-to leading order (N3L0) 3|, In par- 
ticular, the N3L0 model of Ref. [4]] (N3LO-Idaho) reaches the same level of accuracy of the 
CDBonn model. 

The available models for the TNI contain, on the contrary to the two-nucleon interaction 
models, a small number of parameters, usually fixed to reproduce the and/or '^He binding 
energies and, in some cases, the nuclear matter equilibrium density. Among the different 
existing models, we quote only those ones of the Urbana and Tucson-Melbourne families. 
Although constructed within different frameworks, these two families of potentials have 
shown to give similar results, when used in conjunction with a given two-nucleon interaction 
model. Therefore, we have considered the Urbana IX (UIX) TNI in conjunction with both 
the AV18 and N3LO-Idaho two-nucleon interaction models. Finally, it should be noticed 
that within the EFT approach mentioned above, also TNIs appear at the next-to-next-to 
leading order (N2L0) |6|. In particular, we will consider the local version of this N2L0 TNI, 
as given in Ref. 

More recently, a new class of two-nucleon interactions has been obtained {Viow-k po- 
tentials). With the purpose of eliminating from the semi-phenomenological high precision 
two-nucleon potentials the high-momentum parts, the two-nucleon Hilbert space has been 
separated into low- and high-momentum regions and the renormalization group method has 
been used to integrate out the high-momentum components above a cutoff A [8] . The value 
for A is typically fixed in A > 2 systems, for example so that the triton binding energy is 
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reproduced. 

At this point, a crucial issue is to test the model for the nuclear Hamiltonian studying 
A > 4 bound states and A > 3 scattering states. In the present work, we focus our 
attention to the A = 3 scattering problem, which has been the object of a large number of 
investigations . Traditionally, the A = 3 scattering problem with realistic Hamiltonians is 
solved using the Faddeev equations. On the other hand, we have developed in recent years 
a variational approach, based on the expansion of the wave functions on the hyperspherical 
harmonics (HH) basis (for a recent review, see Ref. [l^). This method has proven to 
be very efficient in the description of bound and scattering states in few-nucleon systems. 



In Ref. 



ll| the HH expansion with correlations factors (the correlated HH - CHH - and 



pair-correlated HH - PHH - expansions) has been used to describe A = 3 bound states 



whereas the extension to scattering states has been discussed in Ref. [12|. The inclusion of 
correlation factors was motivated by the short range repulsion of the two-nucleon potential 
which induces particular configurations in the wave function difficult to describe using the 
bare expansion. In fact, in Ref. ISj the HH expansion without correlation factors has been 
used to describe the A = 3 bound state, with the AV18 interaction. The conclusion was that 
a much higher number of states are necessary when the bare expansion is used. The same 
observation has been done in the A = 4 system and is a direct consequence of using local 
interactions, which result to have a strong repulsion at short distances. The implementation 
of the HH method in momentum-space has been done in Ref. 15| for the A = 3,4 bound 
states. This analysis has revealed a much faster convergence of the expansion when non-local 
potentials are considered, even when TNI terms are taken into account. 

The aim of the present work is twofold. First, we want to extend the HH method to 
describe N — d scattering states using either local or non-local potentials. We will show 
that we can apply the method in both configuration and momentum spaces. Second, we 
will present a detailed comparison between the predictions of the different models, local and 
non-local, at low center-of-mass energies, for n — as well as p — scattering. Moreover, 
we will consider the Coulomb potential plu s the magnetic moment (MM) interaction that 



was shown to give sizable contributions 



16| . To our knowledge, this is the first time that 



non-local two- plus three-nucleon potentials are used to describe p — d scattering at very low 
energies. 

The paper is organized as follows: in Sec. [Tll the HH method for the low-energy scattering 
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problem is described, putting more emphasis on those new developments of the method 
necessary in order to use non-local interaction models. In Sec. Illlt the results for the zero- 
energy scattering lengths and low-energy elastic scattering observables are presented and 
discussed. Some concluding remarks are given in Sec. IIV[ 

II. FORMALISM 

In this section we present the HH method for scattering states. The method for bound 

□ 

states has been most recently reviewed in Ref. [10|, and its main characteristics are briefly 
summarized in the following subsection. 

A. The HH Method for Bound States 

The nuclear wave function for the three-body system with total angular momentum J, 
can be written as 

\'^'-'') = Y.^.K'') ' (2-1) 

where is a suitable complete set of states, and /i is an index denoting the set of 

quantum numbers necessary to completely specify the basis elements. 

The coefficients of the expansion can be calculated using the Rayleigh-Ritz variational 
principle, which states that 

(5^^^-^^ = , (2.2) 

where (5c indicates the variation of vl/*^-^^ for arbitrary infinitesimal changes of the linear 
coefficients c^. The problem of determining and the energy E is then reduced to a 
generalized eigenvalue problem, 

(K'' \H-E I m^^^ ) V = . (2.3) 

The main difficulty of the method is to compute the matrix elements of the Hamiltonian H 
with respect to the basis states |\E'^'^^). Usually H is given as a sum of terms (kinetic energy, 
two-body potential, etc.). The calculation of the matrix elements of some parts of H can 
be more conveniently performed in coordinate space, while for other parts it could be easier 
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to work in momentum space. Therefore, it is important that the basis states I^I/^"^^) have 
simple expressions in both spaces. The HH functions indeed have such a property. 

Let us first consider the expression of the HH functions in coordinate space. The internal 
dynamics of a system of three nucleons of identical mass m is conveniently described in 
terms of the Jacobi vectors xip, X2p, constructed from a given particle permutation denoted 
with p, which specifies the particle order i,j, k, and given by 

Xip = ^Jlin - ^{ri + rj)) . (2.4) 
Here p = 1 corresponds to the order 1,2,3. It is convenient to replace the moduli of and 



Xip with the so-called hyperradius and hyperangle, defined as [17 1 



tan0p = ^ . (2.6) 

Note that p does not depend on the particle permutation p. The complete set of hyper- 
spherical coordinates is then given by {p,ni!'^}, with 

= [x,p,X2p;(Pp\ , (2.7) 

and the suffix (p) recalls the use of the coordinate space. 
The expansion states I^E'^"^^) of Eq. (12. ip are then given by 

where fi{p) for / = 1, . . . M is a complete set of hyperradial functions, chosen of the form 



//(P)=7y^^ Lf(7p)e-i^. (2.9) 

Here Lf'\'~fp) are Laguerre polynomials, and the non-linear parameter 7 is variationally 
optimized. As an example, for the N3LO-Idaho potential, it can be chosen in the interval 
6-8 fm-^ 

The functions 3^{G}(fi^^^) are written as 

3 

yiG}in^'^) = Yl ® ® ^]ss] [T2 ^]tt. , (2.10) 

p=i 



where the sum is performed over the three even permutations. The spins (isospins) of 
particle i and j are coupled to 5*2 (T2), which is itself coupled to the spin (isospin) of the 
third particle to give the state with total spin S (isospin T,Tz). The total orbital angular 
momentum L and the total spin S are coupled to the total angular momentum J, J^. The 
functions ¥^^^^^{^11^^), having a definite value of L, L^, are the HH functions, and are written 
as 

YlG^'H'^) = \Ye,{x2p) ® r,,(xip)j N[G] {cos (f)pY' {sin <j)pY' P^^+^'''+^(cos20p) . (2.11) 

Here Y^^{xip) and Y^^{x2p) are spherical harmonics, N^g] is a normalization factor and 
ph+^i2+2 (^gQg20p) is a Jacobi polynomial, n being the degree of the polynomial. The grand 
angular quantum number G is defined as G = 2n + £1 + ^2- The notations [G] and {G} of 
Eqs. (12. lip and (12.101) stand for [(ii^l2]n] and {£1, £2, 5'2, T2, S", T; n}, respectively, and ix 
of Eq. (12. 8p is = {G},1. Note that each set of quantum numbers {£1, £2, -^5 'S'2, T2, 5*, T} 
is called "channel", and the antisymmetrization of 3^{g}(^'''''*) requires £2 + 5'2 + T2 to be 
odd. In addition, £1 + £2 must be even (odd) for positive (negative) parity. To be no- 
ticed that after the sum on the permutation in Eq. (I2.10p . some states inside the sub-space 
spanned by G are linearly dependent. These states have been identified and removed from 
the expansion Q, Q . 

In this work, we have considered modern two-body potential models which act on specific 
spin and angular momentum states of the two-body system. Due to the presence of the sum 
over the permutations in the expression for 3^{g}; a given particle pair is not in a definite 
angular and spin state. However, the HH functions with the grand angular quantum number 
G constructed in terms of a given set of Jacobi vectors Xip, X2p, defined starting from the 
particle order i,j, k, can be always expressed in terms of the HH functions constructed, for 
instance, in terms of a;i(p=i), X2(^p=i). In fact, the following relation holds 

■^[il/2;n](^p'^) = '^^i/2,n;^;,^2'"''^K^4;"'l'-^(P^ ' (2.12) 



where the sum is restricted to the values £[, £'2, and n' such that i[ + £'2 + 2n' = G. The 
coefficients a^^^.p (i relating the two sets of HH functions are known as the Raynal-Revai 



coefficients 18|, and could be computed rather easily using the ortho normality property of 
the HH functions, namely 



Also the spin-isospin states can be recoupled to obtain states where the spin and isospin 
quantum numbers are coupled in a given order of the particles. The result is that the 
antisymmetric functions 3^{g} can be expressed as a superposition of functions constructed 
in terms of a given order of particles i,j, k, each one having the pair i,j in a definite spin 
and angular momentum state. When the two-body potential acts on the pair of particles 
i,j, the effect of the projection is easily taken into account. 

The expansion states of Eq. (12.11) in momentum space can be obtained as follows. Let 
hkip, hk2p be the conjugate Jacobi momenta of the Jacobi vectors, given by 

^kip = y|(pfe - ^(Pi+Pj)) , (2.14) 

Pi being the momentum of the i-th particle. We then define a hypermomentum Q and a set 
of angular-hyperangular variables as 



Q — y^ip + 5 

= [k2p,kp;^p] , (2.15) 

where 

tan (pp = . (2-16) 

k2p 

Then, the momentum-space version of the wave function given in Eq. (12.81) is 

\^i''^'^^) = 9G,i{Q)y{G}i^^'^^) , (2.17) 

where y{G}{^^^'') is the same as 3^{g}(^*'^'*) of Eq. fl2.10l) with Xip —> hip, and 

poo 3 

goAQ) = / dp-^JG+2{Qp) flip) ■ (2.18) 

Jo V 

With the adopted form of fi{p) given in Eq. (12. 9p . the corresponding functions QgAQ) can 



be easily calculated, and they are explicitly given in Ref. 



15|- 



B. The HH Method for Scattering States Below Deuteron Breakup Threshold 

We consider here the extension of the HH technique to describe N — d scattering states 
below deuteron breakup threshold, when both local and non-local interaction models are 
considered. 
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Following Ref. 12], the wave function '^^^^^ describing the N — d scattering state with 
incoming orbital angular momentum L and channel spin 5*, parity vr = (— and total 
angular momentum J, J^, can be written as 

where '^^^'^'^^ describes the system in the region where the particles are close to each other 
and their mutual interactions are strong, while "^^^^^ describes the relative motion between 
the nucleon and the deuteron in the asymptotic region, where the N —d nuclear interaction 
is negligible. The function which has to vanish in the limit of large intercluster 

separations, can be expanded on the HH basis as it has been done in the case of bound 
states. Therefore, applying Eq. fl2.ip . the function vE'^'^'^'^^ can be casted in the form 

|^L5Jj.^^^^^|^jj.^ , (2.20) 

where I^E'^"^"') is defined in Eqs. (12. 8p and (I2.17p in coordinate- and momentum-space, re- 
spectively. 

The function is the appropriate asymptotic solution of the relative N — d 

Schrodinger equation. It is written as a linear combination of the following functions, 

3 

^W = E^W(P) ' (2-21) 

p=i 

where the sum over p has to be done over the three even permutations necessary to anti- 
symmetrize the functions ^\sjJz^ 



1=0,2 



S J J J: 



X [T2 ® ^]tt. . (2.22) 

Here the spin and isospin quantum numbers of particles i and j have been coupled to ^2 and 
T2, with 5*2 = 1, T2 = for the deuteron, Wi{x2p) is the deuteron wave function component 
in the waves Z = 0, 2, yp is the distance between and the center of mass of the deuteron, 
i.e. Up = ^J^Xip, Yi{x2p) and Yiiyp) are the standard spherical harmonic functions, and 
the functions R^ijjp) are the regular (A = R) and irregular (A = /) radial solutions of the 
relative two-body N — d Schrodinger equation without the nuclear interaction. These regular 
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and irregular functions, denoted as J-'iiyp) and Giiyp) respectively, have the form 



Riiy,) ^ g^y,) = {2L + lyAq^+^CMfniyp)^^^^^ , (2.23) 

where q is the modulus of the N — d relative momentum (related to the total kinetic energy 
in the center of mass system by Tcm = 1^) A* being the N — d reduced mass), t] = 2/ie^/g and 
= QVp SI'S the usual Coulomb parameters, and the regular (irregular) Coulomb function 



Fiijj.^p) {Giirj^^p)) and the factor Ciij]) are defined in the standard way jl9|. The factor 
(2L + l)\\q^CL{fi) has been introduced so that JF and Q have a finite limit for g — > 0. The 
function fniyp) = [1 — exp(— 6?/p)]^^+^ has been introduced to regularize Gl at small values 
of yp. The trial parameter h is determined by requiring that fniyp) — > 1 for large values of 
yp, thus not modifying the asymptotic behaviour of the scattering wave function. A value of 
b = 0.25 fm~^ has been found appropriate. The non-Coulomb case of Eq. (12.231) is obtained 
in the limit 0. In this case, -Fl(^7, ^p)/^p and GL{rj,^p)/^p reduce to the regular and 

irregular Riccati-Bessel functions and the factor (2L + 1)!!Cl(?7) ~^ 1 ^ 0. 
With the above definitions, can be written in the form 

^LSJJ. = [^LU^ss'^ls'jj. + Ks.L'S'mus'jj] , (2-24) 



L'S' 

where the parameters Tl-j^g (g) give the relative weight between the regular and irregular 
components of the wave function. They are closely related to the reactance matrix (/C- 
matrix) elements, which can be written as 

fCis,L'S'iQ) = (2^ + 1)!!(2L' + iy.\q'^+'^'^'CL{v)CL'{v)nis,L's'iQ) ■ (2-25) 

By definition of the /C-matrix, its eigenvalues are tan^^^j, 6lsj being the phase shifts. The 
sum over L' and S' in Eq. (I2.24p is over all values compatible with a given J and parity vr. In 
particular, the sum over L' is limited to include either even or odd values since (—1)"^' = vr. 

The matrix elements TZls l'S' i^) linear coefficients occurring in the expansion 

of '^^^^'^'^ of Eq. (12.201) are determined applying the Kohn variational principle [20]] , which 
states that the functional 

[Tiis.L'sM)] = Tiis.L's.(a)-{<':f-\c\K':!l- 

9 



has to be stationary with respect to variations of the trial parameters in ■ Here E is 

the total energy of the system, m is the nucleon mass, and C is chosen so that 

(^LSJjJ-^l^isjjJ ~ (^isjjJ'^I^LSJj^) = 1 • (2.27) 



As described in Ref. |21j, using Eqs. (12.201) and (12.241) . the variation of the diagonal func- 
tional of Eq. (12.261) with respect to the linear parameters leads to the following system 
of linear inhomogeneous equations: 

Y,{K'''\^\K'')^^' = -DIsjjM ■ (2-28) 

Two different terms corresponding to X = R, I are introduced and are defined as 

DlsjjSf^) = {K''\^\^lsjjJ ■ (2-29) 

The matrix elements 'R-is L'S'il) obtained varying the diagonal functional of Eq. (I2.26P 
with respect to them. This leads to the following set of algebraic equations 

T^LS,L" S"il)^ L' S' ,L" s" = Yls,l'S' , (2.30) 

L"S" 

with the coefficients X and Y defined as 

Yls,l'S' = -{^hjj. + ^c'''"''\mvs'j.j) ■ (2.31) 

Here jg the solution of the set of Eq. (12.281) with the corresponding term. A 

second order estimate of 'R-LS,L's'il) given by the quantities ['J^'[s,L's'il)]j obtained by 
substituting in Eq. (12.261) the first order results. Such second-order calculation provides a 
symmetric reactance matrix. This condition is not a priori imposed, and therefore it is a 
useful test of the numerical accuracy. 

In the particular case of g = (zero-energy scattering), the scattering can occur only in 
the channel L = and the observables of interest are the scattering lengths. Within the 
present approach, they can be easily obtained from the relation 

^''^'^aNd = -\iran'^j,ojiq) ■ (2.32) 
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An alternative way to solve the scattering problem, used when g 7^ 0, is to apply the 



2l|. In this way, the Kohn 



complex Kohn variational principle to the iS-matrix, as in Ref. 
variational principle of Eq. (12.261) becomes 

[Sis,L'S'i<l)] = Sis,L'S'{<i) + . (2.33) 

Here 

with ^i^J^^^ given in Eq. ^OTM and 

p=l 

^tsjjM = i^^fsjjSp)-^isjjSp)) 

+ Y.^is,uAl){^^L'S'jjSP) + ^us'jjSP))- (2-35) 

US' 

The functions ^\sjjSp) same as in Eq. (12.221) . with R^{yp) = FL{rj,^p) / and 

^iiVp) = fR{yp)G ^p) /^p- Note that, with the above definition, the reactance /C-matrix 
elements can be related to the iS-matrix elements as 

JCis^L'S'iQ) = H)[S-UL'S'iQ) - Sll'Sss'] [Sis,L'S'iQ) + Sll'Sss']-' ■ (2-36) 

The calculation involving has been performed with the HH expansion in 

coordinate- or in momentum-space, depending on what is more convenient, as it has been 
explained for the bound state in the previous subsection. Some difficulties arise for the cal- 
culation of the potential energy matrix elements which involve ^lsjj^^ ^-^^ {'^'^.'^''I^I^lsjjJ 
present in Eq. (Km . and {^l's'jj. + '^c'^''^'^'''^' I^I^lsjjJ of Eq. (1^]) . with A, A' = R, I. 
In the present work, we consider both two- and three-nucleon interactions, and therefore 

^ = E^^^ + E (2.37) 

i<j i<j<k 

We first focus on the two-body contribution. Due to the antisymmetry of the wave function, 
the following relation holds 

{^\J2v,,\n'LsjjM) = mvul^lsjjM) > (2-38) 

where |$) can be either of Eq. fl2.20p or ^l's'jjSp') ofEq. fl2.22p . depending on which 

term among D'^, X, and Y is considered. The potential Vu acts on the particle pair 1,2 
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with total angular momentum j, and orbital angular momentum and spin quantum numbers 
^12 (o^ the bra) and iu, Si2 (on the ket), and can be written as 



Vl2 = vl2ix2(j,=i),X2(p=iy,i'i2,s'i2, ^12,812) 



(2.39) 



in coordinate-space, and 



V12 = ^l2(^2(p=l)5 ^2(p=l); ^'125 s[2Jl2, S12) 



(2.40) 



in momentum-space, where X2(p=i) and A;2(p=i) are the moduli of the vectors defined in 
Eqs. (12 ■4p and (I2.14p . respectively. When local potential models are considered, such as the 
AV18, then 



(^2{p=l); ^2(p=l); ^'l25 "5'l2; -^12) -312) "^^12(^2 (p=l); ^12) ■5'l2) ^12) ■Sl2)'^(a^2(p=l) ~x'2^p^i^) . 



(2.41) 

The first difficulty that needs to be overcame arises from the fact that when the Vu operator 
acts on ^lsjjSp ■'-)' particle pair 12 does not have a well definite orbital and spin 
angular momenta. However, the following relation holds: 

Wi{x2(p^l)) RLiVp^l) Yi{x2(p^l))®YL{yp^i) = J]i^;L^?,(a;i(p=i),X2(p=i)) 

' ^ l',L' 

X h'(a:2(p=i)) ®>L'(l/p=i)l (2.42) 

where A, are the total orbital angular momentum and its third component. The functions 

^iL^a'(2^i(p=i)>2;2(p=i)) are given by 



^Il'^L' (^l{p=l) > X2(p=l) ] 



dx2{p=i)dxii^p=i) 



Y{:{x2i^p=i)) ®Yl,{yp=i) 



A,A^ 



X wi{x2[p^i)) Rlivp^i) Yi{x2(^p^i))®YL{yp^i) _ . (2.43) 



A,A^ 



Once the functions F;'^'^J^!(a:i(p=i), X2(p=i)) have been calculated and the spin-isospin states 
have been also properly recoupled, the effect of the projection operator in V12 is easily taken 
into account. 

A second difficulty arises in the calculation of the potential matrix element, when non- 
local potentials expressed in momentum-space are used. On the contrary to the core part 
of the scattering wave function '^^^^^'^ , which can be alternatively expressed in coordinate- 
or in momentum-space, the asymptotic states ^^sjj^ ^'^ ^'^^ have an easy expression in 
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momentum-space, and are more conveniently expressed and used in coordinate-space. This 
is especially true when the Coulomb interaction is considered, as for the p—d case. Therefore, 
we have decided to perform the Fourier transform of the potential f i2(^', k; i', s', i, s) to work 
only in coordinate space, namely 

vi,{r',r-i',s',i,s) = Pdk k'^dk'je{k'r')vi^{k',k-i',s',i,s)j,{kr) , (2.44) 

where je{kr) and je{k'r') are the standard spherical Bessel functions. The integrations over 
k and fc', which run from to oo, are easily performed when the potential model considered 
does not have a high-momentum tail, but goes rapidly to zero at rather low values of k and 
k' . This is true for the N3LO-Idaho and Vio-u^ potential models, but not for the CDBonn. 
Since the main goal of the present work is to perform a first test of the applicability of the 
HH method to the A = ?> scattering problem using non-local realistic interactions, only the 
N3LO-Idaho and Viou^ two-body potentials have been considered. 

Some remarks about the calculation of the three-body contribution to the potential energy 
operator of Eq. fl2.37p are in order. The TNIs considered in the present work are the Urbana 
IX jsl (UIX) and the N2L0 [7] potentials. The first one is used in conjunction with both the 
AV18 and N3LO-Idaho two-nucleon interactions. In the second case, the parameter in front 
of the spin-isospin independe nt p art of the UIX TNI has been rescaled by a factor of 0.384 to 



fit the triton binding energy 22| (UIXp). The N2L0 TNI has been used only in conjunction 
with the N3LO-Idaho potential model. All these TNIs are local potentials, and have a well 
defined operatorial structure. Therefore, the projection procedure of Eqs. fl2.42p an d (12.431) 



is not needed and the present approach follows the footsteps of the PHH technique 
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23| 



III. RESULTS 

In this section we present our results for n — d and p — d scattering observables at 
center-of-mass energies below deuteron breakup threshold. The interaction models which 
have been used are the AVIS and the N3LO-Idaho two-nucleon, and the AV18/UIX, N3L0- 
Idaho/UIXp and the N3LO-Idaho/N2LO two- and three-nucleon interactions. Note that 
the AV18 and AV18/UIX results are the same as those ones first obtained in Ref. 2J], using 
the PHH expansion. We have considered also the Viow-k model, obtained from the AV18 
two-nucleon interaction with a cutoff parameter A equal to 2.2 fm~^. The cutoff param- 
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eter has been chosen so that the triton binding energy is 8.477 MeV, when the complete 
electromagnetic interaction is used, including neutron charge distribution and MM interac- 
tion effects. On the other hand, when no electromagnetic effects are considered, the triton 
binding energy has been found to be 8.519 MeV. In the scattering problem, only the point 
Coulomb interaction has been considered, except when differently indicated. 

Before presenting the results for the considered low-energy N — d observables, we discuss 
the pattern of convergence for some representative quantities, i.e. the n — d doublet zero- 
energy scattering length ^a„d and the p — d = 1/2+, 1/2^ phase shifts and mixing angles 
at center-of-mass energy Ecm = 2.0 MeV, calculated with the N3LO-Idaho two-nucleon 
interaction model. The angular momentum-spin-isospin channels considered for J'^ = 1/2"*" 
and l/2~ are given in Tables ID and [ITl respectively. The notation is the same as in Eq. (12.111) . 
To be noticed that the scattering channels in the case of J'^ = 1/2^ are ordered for increasing 
values of ii + £2- This is true also for all the channels here considered, except those for 
J'^ = 1/2+ (see Table [T]), where the ordering respects an "historical choice", first done in the 
case of the three- nucleon bound state in Ref. 25 1. 

In Table IIIII we present the results for ^a„d and p — d = 1/2+, 1/2^ phase shifts and 
mixing angles [Slsj^ ^) at Ecm = 2.0 MeV, for increasing values of the Laguerre polyno- 
mials M in the hyperradial functions (see Eqs. (12.91) and (I2.18P ). All the 23 (25) angular 
momentum-spin-isospin channels of Table HI f ln|) are considered for J'^ = 1/2+ (l/2~), and 
HH functions up to grand angular momentum G = 20 (21) for all the channels have been 
included. From inspection of the table, we can conclude that the use of M = 28 is enough 
to reach an accuracy of at least 0.002 fm for the scattering length and four significant digits 
for the phase shifts and mixing angles. In fact, for other p — d scattering channels at some 
of the considered values of Ecm, even M = 24 and M = 20 has been found enough to reach 
the same degree of accuracy. 



To study the convergence on the HH expansion, as it has been done in Ref. [Ijl, we 
have separated the HH functions into classes having particular properties and we have taken 
into account the fact that the convergence rates of the different classes are rather different. 
For instance, we expect that the contribution of the HH functions with lower values of 
^i2,Q = £i,a + ^2,a to bc the most important. Therefore, for all the J'" scattering states, 
except J'^ = 1/2+, the different classes are classified with increasing value of £12,0, up to 
£i2,a < 6, and among those ones with the same ii2,a, we have included first the contributions 
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TABLE I: Angular momentum, spin and isospin quantum numbers for the first 23 channels 
sidered in the expansion of the J'^ = 1/2+ core wave function. 



a 




ha 


La 






Sa Ta 


1 











1 





1/2 1/2 


2 














1 


1/2 1/2 


3 





2 


2 


1 





3/2 1/2 


4 


2 





2 


1 





3/2 1/2 


5 


2 


2 





1 





1/2 1/2 


6 


2 


2 


2 


1 





3/2 1/2 


7 


2 


2 


1 


1 





1/2 1/2 


8 


2 


2 


1 


1 





3/2 1/2 


9 


1 


1 





1 


1 


1/2 1/2 


10 


1 


1 


1 


1 


1 


1/2 1/2 


11 


1 


1 


1 


1 


1 


3/2 1/2 


12 


1 


1 


2 


1 


1 


3/2 1/2 


13 


1 


1 











1/2 1/2 


14 


1 


1 


1 








1/2 1/2 


15 


2 


2 








1 


1/2 1/2 


16 


2 


2 


1 





1 


1/2 1/2 


17 


3 


1 


2 


1 


1 


3/2 1/2 


18 


1 


3 


2 


1 


1 


3/2 1/2 


19 














1 


1/2 3/2 


20 


1 


1 





1 


1 


1/2 3/2 


21 


1 


1 


1 


1 


1 


1/2 3/2 


22 


1 


1 


1 


1 


1 


3/2 3/2 


23 


1 


1 


2 


1 


1 


3/2 3/2 
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TABLE II: Same as Table HI but for the first 25 channels considered in the expansion of the 
J'" = l/2~ core wave function. 



a 


ha 


ha 


La 


-520 


T2a 


Sa Ta 


1 


1 





1 


1 





1/2 1/2 


2 


1 





1 





1 


1/2 1/2 


3 


1 





1 


1 





3/2 1/2 


4 





1 


1 


1 


1 


1/2 1/2 


5 





1 


1 








1/2 1/2 


6 





1 


1 


1 


1 


3/2 1/2 


7 


2 


1 


1 


1 


1 


1/2 1/2 


8 


2 


1 


1 








1/2 1/2 


9 


2 


1 


1 


1 


1 


3/2 1/2 


10 


2 


1 


2 


1 


1 


3/2 1/2 


11 


1 


2 


1 


1 





1/2 1/2 


12 


1 


2 


1 





1 


1/2 1/2 


13 


1 


2 


1 


1 





3/2 1/2 


14 


1 


2 


2 


1 





3/2 1/2 


15 


3 


2 


1 


1 





1/2 1/2 


16 


3 


2 


1 





1 


1/2 1/2 


17 


3 


2 


1 


1 





3/2 1/2 


18 


3 


2 


2 


1 





3/2 1/2 


19 


1 





1 





1 


1/2 3/2 


20 





1 


1 


1 


1 


1/2 3/2 


21 





1 


1 


1 


1 


3/2 3/2 


22 


2 


1 


1 


1 


1 


1/2 3/2 


23 


2 


1 


1 


1 


1 


3/2 3/2 


24 


2 


1 


2 


1 


1 


3/2 3/2 


25 


1 


2 


1 





1 


1/2 3/2 
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TABLE III: n — d doublet scattering length ^a^d in and p — d J'^ = 1/2^, l/2~ phase shifts 
Slsj and mixing angles e at Ecm = 2.0 MeV, calculated with the HH technique using the N3L0- 
Idaho interaction model, for increasing values of the Laguerre polynomials M. All the channels of 
Tables |T] and |TI] are included with grand angular momentum for all the channels set equal to 20 for 
J"" = 1/2+ and 21 for J"" = 1/2-. 

M = 4M = 8M = 12M = 16M = 20M = 24M = 28 
^Qnd 3.029 1.630 1.329 1.259 1.240 1.234 1.233 
(5n 1 1 -3.611 -3.583 -3.572 -3.570 -3.570 -3.569 -3.569 
5. 3 1 -43.28 -34.69 -32.41 -31.96 -31.82 -31.78 -31.77 

^'2'2 

ei+ 0.525 0.975 1.150 1.189 1.201 1.205 1.206 

2 

5, 1 1 -8.270 -7.756 -7.608 -7.581 -7.576 -7.575 -7.575 

' 2 ' 2 

5, 3 1 20.82 21.73 21.97 22.00 22.00 22.01 22.01 

' 2 ' 2 

ei- 4.947 5.542 5.628 5.641 5.643 5.644 5.644 



of the HH functions with lower £2.0- Finally, the Tq, = 3/2 states are considered. With these 
criteria, in the J'" = l/2~ case, the channels have been classified in 6 classes, including 
channels 1-3, 4-6, 7-10, 11-14, 15-18, and 19-25 of Table UTi respectively. In the case of 
J'^ = 1/2+, the classification follows the footsteps of Ref. 25|, and therefore the channels 
have been classified in 5 classes, including channels 1-3, 4-8, 9-12, 13-18, and 19-23 of 
Table [H respectively. We have then called Gi, for each class i, a number such that each state 
of class i has the grand angular momentum G < Gi, and we have increased Gi until we have 
reached convergence. Then, keeping Gi fixed at this value, we have included the states of the 
following class, and increased Gi+i again until we have reached convergence. The results for 
the zero-energy scattering length and the low-energy phase shifts and mixing angles obtained 
with this procedure are given in Table HVl for J'^ = 1/2+ andlVlfor J'^ = 1/2^. Here, M = 28 
Laguerre polynomials in the expansion of the hyperradial function are included, and again 
the N3LO-Idaho two-nucleon potential is used. 

From the cases presented in the Tables, and as well as for all cases taken in consideration, 
we can observe that (i) the last classes of channels, corresponding to the = 3/2 states, 
give sizable contributions to the p — d phase shifts and mixing angles, but negligible ones 
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to the n — d ones, (ii) The = 1/2 channels with the largest values of £12,0 (fourth class 
for J'^ = 1/2+ and fifth one for J'^ = 1/2^) give negligible contributions. This implies that 
^12,0 < 6 is enough to have accurate results, (iii) The convergence with respect to the grand 
angular momentum for the first class is the most problematic and it depends noticeably 
on the interaction. For example, for J'^ = 1/2+, when the non-local potential N3LO-Idaho 
is used, values of Gi up to 80 have been found to be necessary (see Table HVT) . On the 
other hand, in the case of the local AV18, we have verified that within the HH expansion 
(i.e. without the correlation) Gi = 160 is needed to reach the same degree of accuracy. 
This is related to the fact that the AV18 potential is more repulsive at short interparticle 
distances, and therefore the corresponding wave functions in that region are more difficult to 
be constructed. In fact, when the calculation is performed using the Viom-k potential model, 
which is very soft at short interparticle distances, it is sufficient to set Gi = 40. A completely 
identical pattern of convergence is found for all other J'^ waves, (iv) The convergence of the 
other classes is usually faster than for the first class, as it is evident for the cases reported 
in Tables HV] and |Vl For the J'" = 1/2+ case, we obtain convergence with just G2,3,4,5 = 20. 
For J'^ = 1/2^, we have to consider fairly large values of G only for the fourth class (up 
to G4 = 51), since the channels belonging to this class (the channels 11-14 as reported 
in Table [IT]) are needed to describe pairs in orbital angular momentum £2 = 2. Namely, 
together with the channels of the first class, they are needed to have a good descriptions 
of the pairs in the deuteron waves. We have also found that the convergence rate of these 
classes does not depend much on the non-local interaction model. For example, with the 
Viou^k potential, convergence is achieved with ^2,3,4,5 = 20 for J'^ = 1/2+ and 6*2,4 = 31, 
^3,5,6 = 21 for = 1/2^. However, note that for the AV18 potential model, we need to 
set G2 = 90, ^3,4,5 = 40 for = 1/2+ and G2 = 61, 63,6 = 41, G4 = 91, G5 = 21 for 
J'^ = 1/2^. A similar pattern of convergence has been found for all the calculated quantities. 
From now on, all the results which will be presented have been obtained at convergence in 
the basis expansion. 

The results for the n — d and p — d doublet and quartet scattering lengths are given in 
Table |VT] and are compared with the available experimental data 26|, |27|]. The results for the 
AV18 and AV18/UIX have been taken from Ref. Q- Comparing the theoretical and exper- 
imental results for ^a„(i and ^a„d5 we can conclude that '^ami is very little model-dependent 
(as well as ^ctpd), and there is a satisfactory agreement between theory and experiment. On 
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TABLE IV: Convergence of the n — d doublet scattering length '^Und in and p — d J'" = 1/2"*" 
phase shifts 5l5J and mixing angles e at Ecm = 2.0 MeV, corresponding to the inclusion in the 
wave function of the different classes in which the HH basis has been divided. The N3LO-Idaho 
interaction model is used, and the number of Laguerre polynomials included is M = 28. 



Gi 


G2 




G4 G5 


2 


and 


^0 1 1 

' 2 ' 2 


5n 3 1 
' 2 ' 2 


f 


"1 + 
2 


50 










1 


245 


-3.577 


-32.11 


1 


248 


60 










1 


243 


-3.577 


-32.09 


1 


248 


70 










1 


242 


-3.577 


-32.08 


1 


248 


80 










1 


242 


-3.577 


-32.08 


1 


248 


80 


16 








1 


112 


-3.572 


-31.18 


1 


240 


80 


20 








1 


112 


-3.572 


-31.17 


1 


238 


80 


20 


16 






1 


100 


-3.569 


-31.09 


1 


239 


80 


20 


20 






1 


100 


-3.569 


-31.09 


1 


239 


80 


20 


20 


16 




1 


099 


-3.569 


-31.09 


1 


239 


80 


20 


20 


20 




1 


099 


-3.569 


-31.09 


1 


239 


80 


20 


20 


20 


16 


1 


099 


-3.569 


-31.04 


1 


241 


80 


20 


20 


20 


20 


1 


099 


-3.569 


-31.04 


1 


241 



the contrary, ^a^d is strongly mo del- dependent, and only the inclusion of the TNI brings the 
theoretical value close to the experimental one. However, some disagreement still remains. 



and the recent measurement of Ref. 



27( 1 is not well described by any of the potential models 



considered. Though, the N3LO-Idaho/UIXp and N3LO-Idaho/N2LO models give slightly 
better results. Note that the AV18/UIX results obtained including also MM interaction 
effects are ^a„d = 0.590 fm and ^and = 6.343 fm. Finally, the Vion^ results are in remark- 
able disagreement with the experimental data, and a sizable difference from the AV18/UIX 
results is also observed. Therefore, even when the cutoff parameter of the Viou^ interaction 
model is fixed to reproduce the triton binding energy, the doublet scattering length is not 
well reproduced. This observation seems to suggest that the S-wave sensitive scattering 
observables, like the scattering lengths, are not properly described by simply increasing the 
attraction, but a right balance between attraction and repulsion of the nuclear force has to 



19 



TABLE V: Same as Table HVl but for p — d J"^ = 1/2 phase shifts Sisj <ii^d mixing angles e at 
Erm = 2.0 MeV. 





G2 


G3 


G4 G5 




5^ 1 1 
2 ' 2 


^1, 


3 1 
2 ' 2 


£1- 


61 












-7.416 


21 


24 


5.544 


71 












-7.413 


21 


25 


5.545 


81 












-7.412 


21 


25 


5.545 


91 












-7.411 


21 


25 


5.545 


91 


11 










-7.382 


21 


53 


5.619 


91 


21 










-7.380 


21 


55 


5.622 


91 


31 










-7.379 


21 


55 


5.622 


91 


31 


15 








-7.367 


21 


77 


5.704 


91 


31 


21 








-7.367 


21 


77 


5.705 


91 


31 


21 


31 






-7.372 


22 


02 


5.799 


91 


31 


21 


41 






-7.370 


22 


03 


5.798 


91 


31 


21 


51 






-7.369 


22 


04 


5.798 


91 


31 


21 


51 


15 




-7.369 


22 


04 


5.798 


91 


31 


21 


51 


21 




-7.369 


22 


04 


5.798 


91 


31 


21 


51 


21 


15 


-7.342 


22 


05 


5.818 


91 


31 


21 


51 


21 


21 


-7.340 


22 


05 


5.819 



be reached. Such a balance cannot be achieved with just one parameter, as the cutoff A of 



28|. 



the Viou^k interaction. Further analysis of these aspects is currently underway 

The p — d elastic scattering observables have been studied at different values of center- 
of-mass energy Ecm- Since we have considered several interaction models, we first focus 
our attention on the two-nucleon only models, i.e. the AV18 and the N3LO-Idaho. The 
differential cross section da/dQ, the proton vector analyzing power Ay, the deuteron vector 
and tensor analyzing powers iTu, T20, T21 and T22, as function of the center-of-mass arigle 



9r 



30 



are given 



31 



32 



33, 



in Figs. [H El [31 m [5] and [6l respectively. The data are taken from Refs. [29|, 



34l |. as indicated in detail in the figure captions. By inspection of the figures. 



we can observe that: (i) theory and experiment are in disagreement for the Ay and iTu 
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TABLE VI: n — d and p — d doublet and quartet scattering lengths in fm calculated with the HH 
technique using different Hamiltonian models. 



Interaction 


O-nd 


O-nd 


2„ 4„ 
"pd ""pd 


AV18 


1.275 


6.325 


1.185 13.588 


AV18/UIX 


0.610 


6.323 


-0.035 13.588 


N3LO-Idaho 


1.099 


6.342 


0.876 13.646 


N3LO-Idaho/UIXp 


0.623 


6.343 


-0.007 13.647 


N3LO-Idaho/N2LO 


0.675 


6.342 


0.072 13.647 


^low—k 


0.572 


6.321 


-0.001 13.571 


Exp. [26] 


0.65±0.04 


6.35 ±0.02 




Exp. \27} 


0.645±0.003±0.007 







observables (the well-known "A^-puzzle" 2S, ISSj); (ii) no differences between the AV18 
and the N3LO-Idaho curves can be seen for the differential cross sections; (iii) the N3L0- 
Idaho curves are systematically closer to the data than the AV18 ones for the polarization 
observables, especially for Ay and iTn. The reason of this behaviour is well known 16| and 
is related to the MM interaction. In fact, the AV18 potential model has been constructed 
keeping the electromagnetic interaction separated from the nuclear one. The electromagnetic 
interaction includes the MM one, as well as higher-order corrections to the pp Coulomb 
potential as two-photon exchange, Darwin-Foldy and vacuum polarization terms. The MM 
interaction effects are known to be sizable in N—d elastic scattering 16j. On the contrary, the 
N3LO-Idaho potential model keeps as electromagnetic interaction only the point Coulomb 
potential and MM effects are indirectly included in the nuclear part of the interaction by 
the fitting procedure. From this observation, we can guess that the results obtained with 
the two-nucleon potentials AVIS and N3LO-Idaho should be comparable when the AV18 
calculation includes also the MM effects. To verify this hypothesis, we have calculated the 
p — d elastic scattering observables at two values of Ecm, 1-33 and 2.0 MeV, using the AV18, 
AV18+MM, and N3LO-Idaho potential models. The results are given in Figs. [7] and [H 
respectively. From inspection of the figures, we can notice that the AV18+MM results for 
the Ay and iTu vector polarization observables are larger than the AV18 alone ones in 
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FIG. 1: p-d differential cross section for Ecm = 0.266,0.431,0.666, 1.33, 1.66 and 2.0 MeV calcu- 



lated with the AV18 (dashed lines) and the N3LO-Idaho (solid lines) two-nucleon potential models. 

29( at Ecm = 0.266 MeV, from Ref. [30] at Ecm = 0.431 MeV, from Refs. (29 1 



Data are from Ref. 



(solid circles), [SlJ (empty circles), and 



32] (empty squares) at Ecm = 0.666 MeV, from Refs. 



33( 1 (solid circles), anc 



(empty squares ~ Ep = 1.993 MeV) ,^ 
at Ecm = 1-33 MeV, from Refs. ^] (solid circles) and 



Ecm = 1-66 MeV, from Refs. 
(empty circles) at Ecm = 2.0 MeV. 



33] (solid circles) 



32| 



3J (empty circles - Ep = 2.08 MeV) 



3|] (empty circles - = 2.53 MeV) at 



32] (empty squares - Ep = 2.995 MeV), and [34 1 



the maximum region, and that the AVIS+MM and N3LO-Idaho curves are quite close to 
each other for all the observables considered. Although this analysis should be performed 
systematically at any value of Ecm and for any observable, given the conclusions of Ref. 16| , 
it can be expected that a similar behaviour still holds. Therefore, we can conclude that the 
non-local N3LO-Idaho and the local AV18 two-nucleon interactions give similar results once 
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FIG. 2: Same as Fig. [H but for the proton vector analyzing power Ay. Data are from Ref. [30t] at 
E^^ = 0.431 MeV, from Ref. [31] at Ecm = 0.666 MeV, from Ref. [33] at Ecm = 1-33, 1.66 and 2.0 
MeV. 

the MM effects are included in the AV18 calculation. For this reason, we have chosen to use 
the N3LO-Idaho two-nucleon interaction model in the continuation of our study. 

In order to have a meaningful comparison with the data, the TNI cannot be neglected in 
the calculation. Therefore, we present in Figs. [9l [10], [H] [12], [13] and [TH the results for the dif- 
ferent observables, obtained with the N3LO-Idaho two-nucleon, and the N3LO-Idaho/UIXp 
and N3LO-Idaho/N2LO two- and three-nucleon interaction models. From inspection of the 
figures, we can observe that the TNI effects are sizable, especially for the polarization ob- 
servables, and the N3LO-Idaho/N2LO potential model gives a slightly better description 
of the data than the N3L0- Idaho/ UIXp one. In particular, it is interesting to notice that 
the Ay and iTu observables are better described at every value of Ecm, except for iTu at 
Ecm = 1-66 MeV, although even in this case all the curves are very close to each other. 
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FIG. 3: Same as Fig. [H but for the deuteron vector analyzing power iTn. Data are from Ref. 3C| 



331 ] at Ecm = 1-66 



at Ecm = 0.431 and 1.33 MeV, from Ref. [31] at E^m = 0.666 MeV, from Ref. 
and 2.0 MeV. 

For a better comparison between the different potential models and the data, a analysis 
has been carried only for those observables, except the differential cross section, for which 
the number of data is > 7. In particular, following Ref. jsGj ]. 

xVdatum = -^ ^^-^^^J ^ , (3.1) 

where f^^"^ is the ith datum at center-of-mass angle 9i, Afi is its experimental error, and 
is the theoretical value at the same angle. The results are given in Table IVIII for 
= 0.666, 1.33, 1.66 and 2.0 MeV. The N3LO-Idaho, N3LO-Idaho/UIXp and N3L0- 
Idaho/N2LO interaction models have been considered. From inspection of the table we can 
notice that all the values for x^/datum are comparable, although the ones obtained with the 
N3LO-Idaho two-nucleon interaction are usually higher than the ones obtained with two- 
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FIG. 4: Same as Fig. [H but for t 
at Ecm = 0.431 MeV, from Ref. 
MeV. 



le deuteron tensor analyzing power Ton- Data are from Ref. [3C | 



at Ecm = 0.666 MeV, from Refs. 



331] at Ecm = 1-66 and 2.0 



and three- nucleon interactions, except for the tensor analyzing power T20 and T21. This is 
a well-known and still unclear issue, i.e. T20 and T21 are better described, as the energy 
increases, by two-nucleon only interaction models, even at 30.0 MeV js?]. Among the two- 
plus three-nucleon interaction models, the N3LO-Idaho/N2LO performs slightly better than 
the NSLO-Idaho/UIXp. 

The p — d elastic scattering observables at Ecm = 0.666 and 2.0 MeV have been calculated 
also using the two-nucleon only potential model Viom-k-, obtained from the AVI 8 with a 
cutoff parameter A equal to 2.2 fm~^, as already used for the calculation of the scattering 
lengths. The results are given in Figs. [15] and [161 respectively. Together with the Viou^ 
results, we have shown also the bare AV18 and the AV18/UIX ones. From inspection 
of the figures, we can observe that the Viov^ results are very similar to the AV18/UIX 
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FIG. 5: Same as Fig. [H but for the deuteron tensor analyzing power T21. Data are from the same 
references as in Fig. [H 

ones. This can be understood noticing that the considered observables are sensitive to P- 
and D-wave scattering. The P-wave phase shifts and mixing angles are influenced by the 
UIX TNI attraction term, which is reproduced, within the Viou}~k approach, by fitting the 
cutoff parameter A. In fact, the J'^ = l/2~ phase shifts and mixing angle (^ii i , 5^3 1 , ei-) 

2 2 2 2 2 

obtained at Ecm = 2.0 MeV with the AV18, AV18/UIX and Viou>-k potential models are 
(-7.358,22.11,5.718), (-7.366,22.32,5.835), and (-7.343,22.26,5.811), respectively. From 
this first analysis of Viouh-k results for N — d scattering at low energies, we can conclude that 
the Vion^ and AV18/UIX results are close to each other for observables sensitive to P- and 
D-wave scattering, like vector and tensor analyzing powers. Further work on these aspects 
is currently underway. 

The n — d elastic scattering observables, including differential cross section, neutron vector 
analyzing power Ay, deuteron vector and tensor analyzing powers iTn, T20, T21 and T22, at 
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FIG. 6: Same as Fig. [H but for the deuteron tensor analyzing power T22- Data are from the same 
references as in Fig. [H 

Ecm = 1-33 and 2.0 MeV are given in Fi gs. fTTl and fT8 | respectively. The experimental data 



are from Refs. 



38|,l39 



40| and Refs. |4l|, |42| at E^^ = 1.33 MeV and 2.0 MeV, respectively. 



The different curves are obtained using the N3LO-Idaho, N3LO-Idaho/UIXp and N3L0- 
Idaho/N2LO potential models. From inspection of the figures, we can observe that all the 
curves are very close to each other, especially for the differential cross section da/dQ and the 
tensor analyzing powers T20, T21 and T22, although some small differences are appreciable. 
Moreover, some differences are present for the Ay and iTu vector polarization observables 
at the peak, even if TNI effects are small. Comparing the calculations with the data, we can 
observe that the calculated da/dQ at Ecm = 1-33 MeV is much lower than the measured 
one for large values of the center-of-mass angle 9cm- Such a discrepancy however disappears 
at Ecm = 2.0 MeV. This difference has been observed before and its origin has still to be 
clarified As in thep — ci case, the n — d vector analyzing powers Ay are poorly described 
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FIG. 7: Theoretical results for p — d differential cross section da/dQ, and polarization observables 
Ay, iTii, T20, T21 and T22-, at Ecm = 1-33 MeV are compared to the experimental data. The 
calculation are done using the AV18 (dashed lines), the AV18+MM (dotted-dashed lines), and the 
N3L0- 

MeV), 33(] (solid circles), and |3^ (empty circles - Ep = 2.08 MeV) for the differential cross section, 



daho (solid lines) interactions. The data are from Refs. [34] (empty squares - Ep = 1.993 



and from Refs. 



33] and 



30] for the Ay and iTn polarization observables, respectively. The incident 



proton (deuteron) is Ep = 2.0 MeV {E^ = 4.0 MeV). 

28 





30 60 90 120 150 180 30 60 90 120 150 180 
e [deg] 0^^ [deg] 



0.03 



0.02 



iT 



11 



0.01 - 









0.02 



- T 



20 



- -0.02 



0.04 



30 60 90 120 150 180 30 60 90 120 150 180 
e [deg] e [deg] 





21 





30 60 90 120 150 180 30 60 90 120 150 180 
9 [deg] [deg] 



FIG. 8: Same as Fig. [7] but for E^m = 2.0 MeV. The data are from Refs. 



331 ] (solid circles), 



^1 



(empty squares - Ep = 2.995 MeV), and [34] (empty circles) for the differential cross section, and 



from Ref. 



331 ] for the polarization observables. The incident proton (deuteron) is Ep = 3.0 MeV 



{Ed = 6.0 MeV). 
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FIG. 9: p-d differential cross section for Ecm = 0.266,0.431,0.666, 1.33, 1.66 and 2.0 MeV calcu- 
lated with the N3LO-Idaho (dashed lines), the N3LO-Idaho/UIXp (dotted-dashed lines) and the 
N3LO-Idaho/N2LO (solid lines) two- and three- nucleon interaction models. Data are from the 
same references as in Fig. [TJ 

by the theory in the maximum region, but it should be noticed that the N3LO-Idaho/N2LO 
gives again a better description of the observables than the N3LO-Idaho/UIXp Hamiltonian 
model. 



IV. CONCLUSIONS 



FoUowirig our previous studies on the HH method revisited to work in momentum- 



space id, ll5|], we have implemented our technique to study the N — d elastic scattering 
problem at center-of-mass energies below deuteron breakup threshold, using both local and 
non-local realistic nuclear interactions. Using this method, it is possible to accurately cal- 
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FIG. 10: Same as Fig. [9], but for the proton vector analyzing power Ay. Data are from the same 
references as in Fig. [2l 

culate N — d scattering observables at very low energies, including the contribution from 
the Coulomb potential as well as from higher order electromagnetic terms, such as the MM 
interaction. In particular, it is the first time that nuclear model including non-local two- 
nucleon interactions plus TNIs are used to describe p — d scattering at very low energies. 
We have studied several observables, as scattering lengths, differential cross section, vector 
and tensor analyzing powers, and we have compared our results with the available experi- 
mental data. Our main conclusions can be summarized as follows: (i) the results obtained 
from the local AV18 and the non-local N3LO-Idaho two-nucleon interaction are quite dif- 
ferent from each other, especially for the vector polarization observables Ay and iTu in the 
maximum region, (ii) The differences between AVI 8 and N3LO-Idaho results are strongly 
reduced when the MM effects are included in the AVIS calculation. To be noticed that 
the MM effects are indirectly included in the nuclear N3LO-Idaho interaction, since in the 
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FIG. 11: Same as Fig. [H but for the deuteron vector analyzing power iTu. Data are from the 
same references as in Fig. [3l 

fitting procedure for this model only the point Coulomb interaction between pp is used, 
(iii) Among the TNIs here considered, the N2L0 model performs slightly better than the 
UIX one. The N3LO-Idaho/N2LO results are in fact generally closer to the experimental 
data than the N3L0- Idaho/ UIXp ones, (iv) The Viow^ two-nucleon interaction model has 
also been considered, obtained from AV18 with a cutoff parameter A = 2.2 fm~^, fitted 
to reproduce the triton binding energy. The V^o«^-fc results for those observables sensitive to 
S-wave scattering, such as the scattering lengths, are in strong disagreement with the exper- 
imental data and quite different from the corresponding AV18/UIX ones. On the contrary, 
the results for those observables sensitive to P- and D-wave scattering, such as vector and 
tensor analyzing powers, are very similar to the corresponding AV18/UIX ones. Further 
studies on these aspects are currently underway. We expect to extend the present approach 
to the A = 4 scattering problem below breakup threshold, as already done for zero-energy 
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FIG. 12: Same as Fig. [U but for the deuteron tensor analyzing power T20. Data are from the same 
references as in Fig. [H 
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FIG. 16: Same as Fi g. [151 but at Ecm = 2.0 MeV . The data are from Refs. 



(empty circles) and 



si ,„M Cities. 



^1 



3j] (empty squares) for the differential cross section, and from Ref. [33l] for 



the polarization observables. The incident proton (deuteron) is Ep = 3.0 MeV {Ea = 6.0 MeV). 
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FIG. 17: n — d differential cross section da/dQ, Ay, iTn, T20, T21 and T22, at Ec 



1.33 MeV 



are calculated with the N3LO-Idaho (dashed line), the N3LO-Idaho/UIXp (dotted-dashed line) 



and the N3LO-Idaho/N2LO (solid line) potential models. The experimental data are of Refs. 



(solid circles) and 



^1 



(empty squares - S„ = 2.016 MeV) for da/dVL, and Ref. [40.] for Ay. The 



incident neutron (deuteron) is En = 2.0 MeV {E^ = 4.0 MeV). 



39 



500 



^ 400 





Q I I I I I I I I I I I I I I I I I I I 1^ I I I I I I I I I I I I I I I l\l Q 

30 60 90 120 150 180 30 60 90 120 150 180 



e [deg] 



9 . [deg] 



0.03 



0.02 



iT 



11 



0.01 - 





Q 1^ I I I I I I I I I I I I I I I I M r I I I I I I I I I I I I I I I I I -I Q QO 

30 60 90 120 150 180 30 60 90 120 150 180 ' 



e [deg] 



e [deg] 



21 





30 60 90 120 150 180 30 60 90 120 150 180 
9 [deg] [deg] 



0.04 



FIG. 18: Same as Fig. [T71 but for Ecm = 2.0 MeV. The experimental data are from Ref. 



4l| for 



da/dn, and Ref. [42] for Ay. The incident neutron (deuteron) is En = 3.0 MeV {E^ = 6.0 MeV). 
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